Search for CP Violation in the Decays — )■ K°F° 
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We have searched for CP violation in the decays -> KgP'^ where denotes a neutral 
pseudo-scalar meson that is either a vr", rj, or rj' using KEKB asymmetric-energy e'*'e~ collision data 
corresponding to an integrated luminosity of 791 fb~^ collected with the Belle detector. No evidence 
of significant CP violation is observed. We report the most precise CP asymmetry measurement 

in the decay K'^n° to date: A^""*^^^"" = (-0.28 ± 0.19 ± 0.10)%. We also report the 

first measurements of CP asymmetries in the decays D^' — > A'577 and D'^ — > K%ri': A^p^^^^ = 

(+0.54 ± 0.51 ± 0.16)% and A^"^^^^"' ^ (+0.98 ± 0.67 ± 0.14)%, respectively. 

PACS numbers: 11.30.Er, 13.25.Ft, 14.40. Lb 



The recent evidence for 13° — D'^ mixing [H-Q and the 
corresponding mixing parameters [3] are at the upper 
edge of standard model (SM) predictions @. However, 
large theoretical uncertainties in these predictions limit 
the sensitivity to effects of physics beyond the SM. An al- 
ternative, potentially more promising approach to search 
for new physics (NP) is the study of violation of the com- 
bined Charge-conjugation and Parity symmetries (CP) 
in the decays of charmed mesons [ij . In contrast to mix- 
ing, the expected SM CP violation in the charm sector 
is small 0- 

In this Letter we report time-integrated CP asymme- 
try measurements in the decays 13° — J> KgP'^ 8] where 
P° denotes a neutral pseudo-scalar meson: 7r°, r], or 77'. 
The time-integrated asymmetry, Acp, is defined as 



^CP 



r(D" 



KOPO) 



r(L)0 



r{DO ^ K^P^) + r(L)" K^P^) 



, (1) 



where F is the partial decay width. 

The observed K^P^ final states are mixtures of — > 
K^P^ and D'^ — s- K^P^ decays where the former 
are Cabibbo-favored (CF) and the latter are doubly 
Cabibbo-suppressed (DCS). In the absence of direct CP 
violation in CF and DCS decays, as expected in the SM, 
the CP violation in these processes within the SM is 
generated from mixing and interference of decays with 
and without mixing, which is parameterized by a'"*^ (we 
adopt the symbols used in Ref. [il ) . SM if ° — K'^ mixing 
leads to a small CP asymmetry in final states contain- 
ing a neutral kaon, even if no CP violating phase exists 
in the charm decay. The asymmetry that is expected 



from the SM is measured to be (-0.332 ± 0.006)% 
from if£ semileptonic decays and referred to as A^°p [13] , 



which is reflected in the value of A, 



""^-^ if DCS decay 

contributions are ignored. Since the a'"'^ value expected 
from the SM is at most 0(10"'') [gI, the value of CP 
asymmetry in the decays D'^ KgP'^ within the SM 
is approximately A§p. On the other hand, if NP pro- 
cesses contain additional weak phases other than the one 
in the Kobayashi-Maskawa ansatz interferences be- 



ep 



tween CF and DCS decays could generate 0(1)% direct 
CP asymmetry in the decays D° — )> KgP^ [12]. Physics 
beyond the SM could also induce 0(1)% indirect CP 
asymmetry 0. Thus, observing Acp inconsistent with 
A^p in D° — >• KgP° decays would be strong evidence 
for processes involving physics beyond the SM 0, [l3] ■ 

In addition to Acp measurements, we examine the 
universality of a'"'* in Z)" decays [6| by comparing our 

previous result [4] with the A^p ^ value reported in 
this Letter. Our previously measured^values of direct CP 

and 



violation asymmetries (denoted [g] 
^_ [l3[ are also updated. 



The decay D*+ — > D^T^t is used to identify the flavor 
of the meson from the charge of the low momentum 
pion (referred to as "the soft pion"), 7r+. Thus, we de- 

,D°^K%P' 



termine A 



CP 



by measuring the asymmetry in the 



signal yield 



D' + ^D"7r+ iVrec 



iVrcc 



(2) 
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where iVj-ec is the number of reconstructed decays. The 
measured asymmetry in Eq. ([2]) includes two contribu- 
tions other than Acp- One is the forward-backward 
asymmetry (Apg) due to 7* — Z° interference in e+e^ — > 
cc and the other is a detection efficiency asymmetry be- 
tween positively and negatively charged soft pions {Ae" ) . 
Since we reconstruct the Kg with tt+tt" combinations 
and P*^ with the 77 or 777r^7r~ final states, asymmetries 
in Kg and detection cancel out. Equation ([2]) then 
can be simplified to give 



4^ 



A 



CP 



(3) 



by neglecting the terms involving the product of asym- 
metries, where Acp is independent of all kinematic vari- 
ables, Apg is an odd function of the cosine of the po- 
lar angle of in the center-of-mass system (CMS), 

and Ae ° depends on transverse momentum and polar 
angle of 7r+ in the laboratory frame, while it is uni- 

form in azimuthal angle. To correct for Ae ° we use 
the decays K~tt~^ (referred to as untagged) and 

_D*+ — s> D^Trf K^TT+nf (referred to as tagged), and 
assumes the same ApB for 1)*+ and D'~' mesons. By sub- 
tracting the measured asymmetries in these two decay 
modes, A'^^^'^ss'^'^ and Al^ss<=d^ directly measure the 



Ae" correction factor [1 
rected for A^" (denoted A, 



ij. With a;. 



rcc.corr 



below). 
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we extract Acp and App using 
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The data used in this analysis were recorded at or near 
the T(4S') resonance with the Belle detector [l^ at the 
e+e" asymmetric-energy collider KEKB [3]. The sam- 
ple corresponds to an integrated luminosity of 791 fb~^. 

We apply the same charged track selection criteria that 
were used in Ref. [17[. For soft pions we do not require 



associated hits in the silicon vertex detector, either in 
the z or radial directions 18|. Charged kaons and pions 
are identified by requiring the ratio of particle identifi- 
cation likelihoods [17| to be greater or less than 0.6, re- 
spectively. Kg candidates are reconstructed from pairs 
of oppositely charged tracks that have an invariant mass 



within ±9 MeV/c^ of the nominal K^ mass [i,[l3|- Can- 
didate tt'^ and 77 mesons are reconstructed from 77 pairs 
where the minimum energy of each 7 is required to be 60 
MeV for the barrel and 100 MeV for the forward region of 
the calorimeter [l^l . We require the 77 invariant mass to 
be between 0.11 and 0.16 GeV/c^ for candidates and 
between 0.50 and 0.58 GeV/c^ for 77 candidates. The 
momentum of 77 pairs is required to be greater than 0.5 
GeV/c for both 7r° and rj selections. In order to remove a 
significant tt" photon background contribution under the 
77 signal peak, we combine individual 7 candidates from 
77 — 77 with any other detected 7 in the event. If the 
77 pair invariant mass is in the 7r° mass window, the 7 is 
rejected. Further reduction of the tt^ contribution under 
the 77 signal is achieved by requiring the energy balance 
of the 77 in the 77 decay to be less than 0.8, where the 
energy balance is the ratio of the difference and the sum 
of two 7 energies. Candidate 77' mesons are reconstructed 
in the r]iT^TT~ decay channel. To improve the 77' mass res- 
olution, the four-momentum of the rj is recalculated with 
a nominal rj mass Q constraint. The same minimum 7 
energy requirement used for the D*^ — >■ K^r] selection is 
imposed in the 77' reconstruction. The tt'^ veto, however, 
is not applied since it is found to be unnecessary once the 
invariant mass of the r]'K~^T:~ candidates is required to be 
between 0.945 and 0.970 GeV/c^ and the D'^ mass selec- 
tion requirement, which is described below, is applied. 

The four-momentum of the P" is recalculated from a 
kinematic fit to its nominal mass {9| and combined with 
a Kg to form a _D" candidate. D*^ candidates are recon- 
structed using a soft pion and a D'^ candidate with mass 
in the [1.75,1.95] GeV/c^ (K^tt"), [1.82,1.90] GeV/c^ 
(if^r/), or [1.84,1.89] GeV/c^ [K^r]') interval which de- 
pends on the mass resolution. To remove D*~^ mesons 
produced in B decays, the D*^ momentum in the CMS 
is required to be greater than 2.5 GeV/c. AU selections 
are chosen to maximize iVg/cr^Vg and to minimize the 
peaking backgrounds, where Ns is the signal yield from 
the fit and a^s is the uncertainty in Ns- After applying 
all of the selections described above, the D'^ 
contribution to 13° — -^s""" and the D 
bution to — > ^^577 are found to be negligible in sim- 
ulation studies. Figure [T] shows data distributions of the 
mass difference M{D*) — M{D) for all the decay modes. 

All mass difference signals are parameterized as a sum 
of a Gaussian and a bifurcated Gaussian distributions 
with a common mean. The background is parameterized 
by the form {x — 7Tj^+)"e~^'^^~™T+-', where a and/3 are 
free parameters, r?z^+ is the charged pion mass [9| and 
X is the mass difference. The asymmetry and the sum 
of the _D*+ and D*^ yields are directly obtained from a 
simultaneous fit to the D*"*" and D*~ candidate distribu- 
tions. The common parameters in the simultaneous fit 
are the mean of the Gaussian, the widths of the Gaussian 
and the bifurcated Gaussian, and the ratio of the Gaus- 
sian and the bifurcated Gaussian amplitudes, which are 



X^tt" contri- 
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the same for the M{D*) - M{D) distributions in differ- 
ent KgP° final states and in the shghtly different phase 
spaces of the individual KgP^ modes. Table U lists the 
results of the fits. 
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FIG. 1: Distributions of the mass difference M{D*) - M{D) 
for the studied decay modes. Left plots show the mass differ- 
ence between D*"*" and D" and right plots show that between 
D*~ and I)". Top plots are for the KgTv", middle plots for the 
KgT], and bottom plots for the Kg-rj' final states. Points with 
error bars are the data and the histograms show the results 
of the parameterizations of the data. Hatched areas are the 
background contributions. 



TABLE L The sum (Ns) and the asymmetry (Arec in Eq. ([2}) 
of D*"*" and _D*~ yields from the fits. The uncertainties are 
statistical only. 









Ns 


A,ec (%) 


D*+ 


D\t - 




326303 ± 679 


+0.19 ±0.19 








45831 ± 283 


+1.00 ±0.51 


D'+ 


- 




26899 ± 211 


+1.47 + 0.67 



In order to obtain A^" we first extract ^""(*^sgcd ^gjj^g 
simultaneous fits analogous to those used for the signal 
modes, but instead of the M{D*) — M{D) distribution we 
fit to the M{D) distribution using a similar parameteri- 
zation. The values of ^""j.^'^sged g^j.g evaluated in bins of 
transverse momentum (pl^^o) and polar angle (cos 6*^0) 
of untagged D° — > K^n^ candidates in the laboratory 
frame. The pt and polar angle variables are only weakly 
correlated. Each tagged D* DtTs — >■ Kmrs candi- 
date is then weighted with 1 - 74,"^'J,*'^sgcd ^*+ 
1 + for D*- . Details of the weighting proce- 

dure are described in Ref. After this the remaining 



asymmetry in the tagged decay sample is A^" , which 
is obtained from the simultaneous fits to the weighted 
M{D*) — M{D) distributions with the same parameter- 
ization used in the signal modes, now for bins of P^^+ 

and cos6'''i'. 

■7r+ 

The dominant sources of uncertainty in the Ac " deter- 
mination are the statistical uncertainties in the untagged 
and tagged samples. These are found to be 0.04% and 
0.07%, respectively. Other sources of systematic uncer- 
tainties are found to be negligible. Thus, we assign a 

total systematic uncertainty of 0.08% to the A^' deter- 
mination, obtained by adding the two contributions in 
quadrature. 

The data samples shown in Fig. [T] are divided into bins 

of pJf^ and cos 6''^'' of the 7r+. The A^' correction is 

applied by weighting each D*+ event with 1 — A^^ and 

each D* event with \ + A^" . The weighted mass dif- 
ference distributions in bins of the D*^ polar angle in 
the CMS are fitted simultaneously to obtain the cor- 
rected asymmetry. We fit for the linear component in 
cosO'^f to determine ApB while the Acp component 



is uniform in cos e'^f. Figure [2] shows A 



CP 



and 




< 



|cose^.'!"'| 



|cos6§."-| 



Q O 
< 



< 



°"|cose§.^S| 



|cose^."S| 



Q O 
< 



|cose^.-1 



|COS6™^| 



FIG. 2: Measured Acp (left) and Afb (right) values as a 
function of | cosS^^j'l. Top plots are for A'^tt", middle plots 
for K'gri, and bottom plots for K^rj' final states. The dashed 
curves show the leading-order prediction for A'^g. 



Ap'g as a function of | cos I 



From a weighted aver- 

K°7r° 



age over the |cos6'^'^^| bins, we obtain A^p 

(+0.54 ± 0.51)% 



(-0.28 ± 0.19)%, App^^^" 



and 



A 



(+0.98 ± 0.67)% where the uncertainties 
are statistical only. The x^/d.o.i with respect to the 



CP 
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average over the Icose™^! bins is 5.1/4 {K^tt°), 3.0/3 
(Kgry), or 5.3/3 (K^rj'). The observed Aps values de- 
crease with cos 6*^^^^ as expected from the leading-order 
prediction [20| . The observed deviations from the predic- 
tion are expected due to higher order corrections. Similar 
Afb's were found in previous measurements 13, 14, 2l|. 
The results are validated with toy pseudo-experiments 
and full detector simulation Monte Carlo events. We 
found no systematic deviations from the input values. 

We consider several sources of systematic uncertainty. 
The uncertainty due to the limited size of the tagged and 
untagged samples was discussed above. To estimate the 
systematic uncertainty due to the choice of fitting method 
and parameters, we vary the histogram binnings, fitting 
intervals, and signal and background parameterizations. 
We also consider the systematic uncertainties due to the 
choice of cos 9™+ binning. Finally, we include possible 
effects due to the differences in interactions of K'^ and K''^ 
mesons with the material of the detector as explained in 
Ref. [21|], and assign a systematic uncertainty of 0.06% 
due to this effect. Table HIl summarizes the components 
of the systematic uncertainties. The larger uncertainties 



CP 



due to the choice of fitting method are a 



in A, 

consequence of smaller statistics of these samples. 



TABLE II: Summary of systematic uncertainties in Acp- 



Source 


Kl^" (%) 


K'sV (%) 


kW (%) 


A^ " determination 
Fitting 
cosO^ff binning 
K°/K"-ma,tenal effects 


0.08 
0.02 
<0.01 
0.06 


0.08 
0.12 
0.01 
0.06 


0.08 
0.10 
0.03 
0.06 


Total 


0.10 


0.16 


0.14 



From the total uncertainties shown in Table HIl we ob- 



tain A 



CP - v-0.28 ± 0.19 ± 0.10)%, A^p^^^'' = 

(-h0.54± 0.51 ±0.16)% and A^p^^^''' = (-hO.98 ± 0.67 ± 
0.14)% where the first uncertainties are statistical and 
the second are systematic. Table IIIII summarizes the re- 
sults, current world average @, and A^p. 



we test the universality of a assuming negligible new 
CP violating effects in D° decays to the K^g'n^ final 
state as discussed in Ref. 0). By subtracting A^p from 

^CP^^^''°, we obtain a'"'^ = (+0.05 ± 0.19 ± 0.10)%, 
which is consistent with -^r = (-0.01 ± 0.30 ± 0.15)% 
obtained in Ref. 2]. This is the first experimental test 
of a'"'^ in decays with a sensitivity near 0.3%. By 
averaging the two independent values we obtain a'"*^ — 
(-1-0.03 ± 0.18)%, where the uncertainty includes the sta- 
tistical and systematic errors, and represents the most 
precise value of a™'^ from a single-experiment currently. 
Using the average a'"*^, we also update the values of 

o-DO^K+K- ^^'^ "do^tt+tt- Ref. Ilsl . which are 

(-0.46 ±0.37)% and (+0.40 ± 0.56)% ^respectively. 
The errors include all the uncertainties of input measure- 
ments. 

In summary, we report a search for CP violation in 
the decays D° K'gP^ using a data sample with an in- 
tegrated luminosity of 791 fb~^ collected with the Belle 
detector. We observe no evidence for CP violation. The 
measurement in the decay — K^tt^ is the most pre- 
cise measurement of any CP asymmetry in the charmed 
particle sector to date. We also report the first measure- 
ments of CP asymmetries in the decays Kgi] and 
D'^ — > Kgi]'. Our results are consistent with the SM and 
can be used to place the most stringent constraints on NP 
models arising from the measurements of CP violation in 
the charm sector at present. 
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TABLE III: Summary of the Acp measurements. The first 
uncertainties in the second column are statistical and the sec- 
ond are systematic. The third column shows the world av- 
erage of Acp and the fourth A§p. A^p^^^^ are the first 
measurements, hence no world average of Acp is given in the 
third column. 

Belle (%) Ref. [9] (%) A§p (%) 
A^p^^s'^" -0.28+0.19+0.10 +0.1+1.3 -0.332+0.006 
^D«-^K°v ^_o.54±o.51+0.16 — -0.332+0.006 
^DO_^A'g„' 4^0.98+0.67+0.14 — -0.332+0.006 



Besides the Acp measurements listed in Table IIIII 
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